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ONE-DIMENSIONAL  MODEL  FOR  MUD  FLOWS1 
David  R.  Schamber2,  A.M.  ASCE  and  Robert  C.  MacArthur^,  A.M.  ASCE 


ABSTRACT 

In  this  paper  V  transient,  one-dimensional  model  for  dynamic  flood 
routing  of  mud  flo»«  Is  presented.  The  governing  equations  of  mass  and 
momentum  conservation  Incorporate  laminar  flow  resistance  effects  and 
utilize  a  power  law  expression  t.o  represent  the  cross-sectional  geometry 
of  the  channel.  The  equations  are  solved  by  the  method  of 
characteristics  on  fixed  time  lines  and  program  execution  is  performed 
on  a  micro-computer.  Numerical  results  are  compared  with  published 
experimental  data  for  a  laminar  flow,  dambreak  problem  of  a  viscous  oil. 


INTRODUCTION 

During  the  spring  of  1983,  widespread  landslides  and  debris  flows 
caused  an  estimated  250  million  dollars  in  damage  In  the  state  of  Utah. 
Along  a  thirty-mile  length  of  the  Wasatch  Front  Mountains,  over  ninety 
significant  landslides  and  debris  flows  sent  torrents  of  mud,  debris  and 
water  down  steep  canyons  onto  residential  areas  located  on  alluvial  fans 
at  the  base  of  the  mountains. 

The  ability  to  model  these  types  of  events  is  clearly  needed  and 
will  be  useful  in  preparing  maps  which  delineate  potential  flood  damage 
areas.  The  purpose  of  this  paper  is  to  present  a  one-dimensional 
mathematical  model  which  can  be  used  to  route  a  mudflow  down  a  confining 
channel.  Equations  of  mass  and  momentum  conservation  are  presented, 
with  frictional  resistance  terms,  which  account  for  the  laminar  flow  of 
a  Bingham  plastic  fluid.  The  equations  are  solved  by  the  method  of 
characteristics  on  fixed  time  lines.  To  verify  the  model,  comparison  is 
made  with  experimental  results  of  a  laminar  flow  dambreak  problem. 


GOVERNING  EQUATIONS 

The  flow  is  governed  by  the  equations  of  mass  and  momentum 
conservation  which  are  given  respectively  by  [6] 
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In  wMeh  x  *  coordinate  along  the  channel;  t  *  time;  A  =  cross-sectional 
area  of  flow;  V  *  average  velocity;  B  *  channel  top  width;  y  *  flow 

depth;  =  rate  of  change  of  area  with  x  for  a  constant  depth 

(nonprl static  term);  g  «  gravitational  constant;  S0  a  slope  of  the 
channel  bottom;  and  Sf  *  resistance  slope. 

In  most  hydraulic  applications,  the  flow  Is  turbulent  and  Sf  is 
generally  given  by  Manning's  equation.  The  flow  of  mud  presents  an 
entirely  different  situation.  DeLeon  and  Jeppson  [1]  summarize  the  data 
from  a  number  of  debris  flows,  mud  flows  and  pipe  sludge  flows  and 

conclude  that  the  flow  Is  usually  laminar.  By  fitting  a  line  through  a 

number  of  data  points,  these  authors  postulate  a  power  law  relation 
between  the  Chezy  coefficient  and  the  flow  Reynolds  number.  Jeyapalan 
et  al_.  [2],  in  their  analysis  of  mine  tailing  dam  failures,  develop  an 
expression  for  Sf  by  analyzing  the  laminar  flow  of  materials  with 
Bingham  plastic  fluid  characteristics.  Other  researchers,  [1,5]  have 
noted  a  similar  behavior  for  mud  flows,  which  often  exhibit  plug  like 
flow  with  a  critical  yield  stress.  In  this  work,  the  resistance  term 
for  a  Bingham  plastic  fluid  Is  adopted  [2],  Mathematically, 
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in  which  np  *  plastic  viscosity;  y  *  unit  weight  of  the  fluid;  xy  = 
yield  stress  of  the  fluid;  h  »  hydraulic  depth;  and  R  *  hydraulic 
radius.  The  first  term  on  the  right  hand  side  of  Eq.  3  is  similar  In 
form  to  the  expression  postulated  by  DeLeon  and  Jeppson  [1], 

Equations  1  and  2  are  hyperbolic  In  nature  and  have  the  property 
that,  through  linear  combination,  they  can  be  reduced  to  equations 
Involving  differentiation  In  one  less  direction  than  the  original 
equations  [6],  This  characteristic  form  for  Eqs.  1  and  2  is  given  by 
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in  which  c  *  celerity  of  an  elementary  gravity  wave  is  given  by 
c  -  (f)15  <6> 
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curves  In  Fig.  1,  l.e.,  L-P  and  R-P  are  approximated  by  parabolas  In  the 
x-t  plane.  With  this  approximation,  the  finite  difference  form  of  Eqs. 
5  Is  given  by 
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In  which  A^  *  Ap  =  Ar  *  %  .  The  forward  and  backward  version  of  Eqs.  4 
are  also  written  in  finite  difference  form.  Mathematically, 


(Vp  +  wp)  -  (VL  +  u)L) 
fit 
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in  which 

F*  *  9<S0  -  sf)  ;  ^  VA*  ?  (»  i  c)  %  £  ff  dA  (16) 


The  set  of  four  nonlinear  equations,  Eqs.  12-15,  determines  the 

locations  of  points  L  and  R  as  well  as  Vp  and  yp.  The  variation  of  yj_, 

VL,  yp  and  VR  Is  determined  by  parabolic  interpolation  along  time  line 
tj\  A  simple  search  procedure  assures  that  the  interpolation  nodes 
(^k-l*  xk»  xk+l)  always  straddle  the  points  in  question,  so  that 
extrapolation  is  avoided. 

Eqs.  12-15  are  solved  iteratively  by  Newton's  method  [4].  A  first 
guess  to  the  solution  is  found  by  solving  a  linear  version  of  Eqs.  12-15 
in  which  A.  *  AR  =  1  and  Ap  =  0.  The  equations  are  solved  at  a  number 
of  points  between  xi  and  xn  to  define  the  wave  profile.  At  the 

boundaries  of  the  flow  domain,  if  only  the  velocity  or  depth  is 
specified,  the  remaining  unknown  Is  determined  by  application  of  the 
appropriate  backward  or  forward  characteristic  equation.  For  the  case 
of  advance  on  a  dry  bed, Whltham's  assumption  Is  used,  i.e.,  Vn  *  Vn_i. 

During  the  early  stages  of  flooding,  the  effects  of  boundary 
roughness  and  channel  slope  are  small.  A  solution  which  ignores 

friction  and  slope  is  therefore  used  as  the  initial  condition  from  which 
to  start  the  numerical  solution. 


RESULTS 

The  model  is  compared  with  several  dambreak  experiments  performed 
by  Jeyapalan  et  al.  [3].  In  these  experiments,  oil  is  used  to  simulate 
a  laminar  flow  o7“ a  viscous  fluid.  The  experiments  are  conducted  in  a  6 
foot  long  glass  flume  which  has  a  constant  width  of  1  foot.  The  dam  is 
located  4  feet  from  the  downstream  edge  of  the  flume  giving  a  reservoir 
length  of  2  feet.  Table  I  gives  the  parameters  characterizing  the 


Table  I  Flume  Test  Parameters 


Test  No. 


8 (degrees) 

Y(lb/ft3) 

0 

56 

0 

56 

0 

56 

0.078 

0.078 

0.156 


examples  presented  herein.  The  test  numbers  listed  in  Table  I 
correspond  to  several  of  the  flood  examples  presented  in  [3].  In  Table 
I,  *  depth  of  oil  immediately  behind  the  dam  before  failure;  6  = 
bottom  slope  of  the  flume;  and  for  all  cases  Ty  *  o. 

Results  of  the  numerical  simulation  are  presented  in  Figs.  2-3  and 
compared  with  the  available  experimental  data.  The  agreement  between 
theory  and  experiment  is  generally  good.  The  numerical  algorithm  is 
programmed  in  Fortran  and  executed  on  an  Apple  Macintosh  micro-computer. 
Computation  times  are  on  the  order  of  0.25-0.34  seconds  per  computational 
node. 
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